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ABSTRACT 


DLTS and EPR measurements are reported on aluminum-doped silicon that 
has been irradiated at room temperature with high-energy electrons. Compari- 
sons are made to comparable experiments on boron-doped silicon. Many of the 
same defects observed in boron-doped silicon are also observed in aluminum- 
doped silicon, but we have observed several not observed, including the alu- 
minum interstitial and aluminum-associated defects. Damage production model- 
ing, including the dependence on aluminum concentration, will be presented. 


INTRODUCTION 


The radiation- induced defects in silicon responsible for the deteriora- 
tion of solar cells are mainly impurity-associated defects such as the vacancy- 
oxygen-carbon complex (the K-center) (ref. 1). One way to improve the radia- 
tion tolerance of silicon cells is to remove the undesired impurities, espec- 
ially oxygen and carbon from the starting material and to prevent their intro- 
duction during the cell fabrication processes. Alternatively, a way to re- 
duce the importance of such impurities is to introduce sinks which would act 
as recombination centers for vacancies and self-interstitials to reduce the 
chance of the primary defects interacting with the impurities, or perhaps, to 
introduce an impurity which would trap a vacancy or self-interstitial to form 
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an electronically inactive complex which would not then reduce the minority 
carrier lifetime in the device. So that we may produce such radiation re- 
sistant silicon it is important to understand the damage production mechan- 
isms and defect interactions in electron irradiated silicon. Sutdies which 
lead to defect characterization and identification and an understanding of 
defect interactions, particularly of interstitial defects, are in progress. 
This paper describes some preliminary results of our work in aluminum-doped 
silicon. This material is interesting because unlike the boron interstitial 
which is mobile at room temperature, the aluminum interstitial is not mobile 
at room temperature, and hence, its role in damage production should be dif- 
ferent. 


EXPERIMENTAL RESULTS 


Both EPR and DLTS measurements have been made on electron irradiated 
aluminum-doped silicon to characterize and identify the defects produced. 

In all samples oxygen and carbon impurities were present in concentrations 
of 10 16 -10 18 cm“ 3 . Diodes were fabricated from float zone and Czochralski 
grown silicon of 1 ftcm and 10 Qcm resistivity irradiated with 1.0 and 1.5 
MeV electrons. Most of the diodes were Schottky barrier diodes in which 
only majority carrier traps are observed although we did observe the minority 
carrier traps in diffused diodes made from 1 ftcm material. Figure 1 shows 
the DLTS spectrum for a diffused diode. The spectrum shows the divacancy 
peak (Hg) and the K-center peak (H^ and H 5 ) both of which are also present in 
boron doped silicon. In addition, a peak at E v 4- 0.17 eV (H]_) is observed 
and also a minority carrier peak at E c - 0.29 eV (E^) . The peak Ex occurs at 
nearly the same energy as a boron defect, tentatively identified as Bj 4* 0 ^ 
(ref. 2), but the peak intensity is considerably smaller. This could be an 
aluminum complex and further study is needed to identify this defect. 

Figure 2 shows a DLTS spectrum for 10 ficm aluminum-doped silicon, in this 
case a Schottky barrier diode. Again, we observe but here the intensity 
compared to the divacancy peak is much smaller, indicating that the concen- 
tration of the trap has decreased as has the aluminum concentration. This 
data plus the absence of this defect in boron-doped material leads to the 
conclusion that this is an aluminum associated defect. H 3 is the carbon in- 
terstitial which anneals at room temperature to form the carbon-oxygen- 
vacancy complex as was observed in boron-doped silicon (ref. 2,3). 

Figure 3 shows the DLTS spectrum for a 1 ftcm aluminum-doped Schottky 
diode showing that the peak labeled H 4 and H 5 is resolved into two peaks if 
the pulse width is reduced sufficiently. At low fluence we believe only the 
carbon-oxygen-vacancy complex at E v 4- 0.33 eV (H^) is present, but that a 
defect at E v 4* 0.44 eV (H 5 ) appears as the fluence is increased. This is ob- 
served as a shift of the apparent energy level of the peak toward higher ener- 
gy as the fluence is increased. We believe this second defect is also present 
in boron-doped material where we have measured two different capture cross 
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sections for this peak. We have not yet identified the defect corresponding 
to the H 5 peak. This spectrum also shows a defect at E v + 0.55 eV (Hg) which 
occurs in many aluminum-doped samples. We did not observe this defect in 
boron doped samples which suggests that it too is an aluminum associated de- 
fect. 


Isochronal annealing for 20 minute periods was performed on Schottky 
barrier diodes to study the annealing behavior of these defects. A new DLTS 
peak appeared at E y + 0.48 eV (Hy) which appears after annealing at 100 P C, 
continues to grow after the 145°C anneal and then decreases after the 200°C 
anneal and is essentially gone after a 250°C anneal. Another additional peak 
is observed after annealing at 400°C at E v + 0.22 eV (Hg) which increases 
after a 445 °C anneal after which another very small peak at E v + 0.20 eV also 
appears. We are continuing to study these new peaks. 

EPR measurements in aluminum-doped silicon show that, in addition to the 
divacancy (G 6 ) (ref. 4), and the carbon interstitial (G12) (ref. 5) which 
anneals at room temperature to form the carbon-oxygen-vacancy complex (G15) 
(ref. 3,6), aluminum interstitials (G18) (ref. 7,8) are observed as well as 
aluminum-vacancy pairs (G9) (ref. 8 ). The aluminum interstitial was observed 
to anneal at 200°C by thermal processes forming an aluminum substitutional- 
aluminum interstitial pair (G19,G20) and a third spectrum (G21) (ref. 9). 

The aluminum interstitial anneals by an ionization enhanced process at 100°C 
(Troxell, Chatterjee, Watkins and Kimerling, to be published). 

We have not always observed the A defect in our DLTS samples. Al- 
though its energy position in the band gap is E v + 0.17 eV, the DLTS peak is 
observed at the apparent energy position of E v + 0.26 eV because of the 
strong temperature dependence of the capture cross section (ref. 10). We be- 
lieve it is buried in our carbon-oxygen-vacancy peak and would be more easily 
observed in samples with low carbon and oxygen content. Table I lists the 
defects observed by DLTS measurements in both boron-doped and aluminum-doped 
electron irradiated silicon. 


DEFECT PRODUCTION MODELING 


In a previous paper we presented results of a computer model for defect 
production by electron irradiation in boron doped silicon containing carbon 
and oxygen as impurities and showed that the presence of a defect recombina- 
tion center could reduce damage produced at a given fluence (ref. 11). The 
model used is a simplified one described by means of a system of rate equa- 
tions for defect interactions. Free mobile single vacancies, free mobile 
self-interstitial atoms and immobile divacancies are created at constant 
rates during irradiation. We have calculated defect production in boron- 
doped silicon of different resistivities assuming no other impurities are 
present, and also for aluminum-doped silicon with the same assumptions. The 
primary difference in these two systems is that the boron interstitial, 


187 



created by the Watkins replacement mechanism, is mobile at room temperature, 
whereas the aluminum interstitial, created in the same way, is not. 

We assume the following defect reactions in boron-doped silicon: 


I + V -* recombination 

(1) 

1 + B s - B I 

(ref. 

12) 

(2) 

Bj + V + Bg 



(3) 

B i + E s -> B i + 

B s E 

B I B S 

(4) 

B I + V 2 V + 

B s 


(5) 

V + B s - VE S 



(6) 

1 + VB s - B s 



(7) 

1 + V 2 - V 



(8) 

V + V V 2 



(9) 


where V = single vacancy, I = self-interstitial, Bg = substitutional boron, 

Bj = interstitial boron, and V 2 = divacancy. We assume similar reactions in 
the aluminum-doped silicon where A£g = substitutional aluminum and A^x = 
interstitial aluminum replace Bg and Bj in the equations. We have calculated 
the damage under different irradiation conditions, one with continuous con- 
stant flux and the other with pulsed flux, as in the case of a rotating satel- 
lite. Table II lists the rate equations and Table III lists the numerical 
values of parameters used for the calculations . 


We present no definitive results here since we are still testing and im- 
proving the program. Preliminary results show that at higher doping levels 
the damage component due to V and V 2 at a particular fluence decreases, while 
that due to defects involving the dopant increases as more secondary and 
tertiary defects are produced. In the boron-doped material the major defects 
in the more heavily doped material are the vacancy-boron substitutional pair 
and the boron interstitial-boron substitutional pair, whereas in the aluminum- 
doped material they are the vacancy-aluminum substitutional pair and the alum- 
inum interstitial. This was expected since the aluminum interstitial is not 
mobile at room temperature and cannot migrate to an aluminum substitutional 
site. The calculations show equal amounts of damage in both types of mater- 
ials however. Figures 4 and 5 show calculations for two different aluminum 
concentrations. The effect of each of these defects on the minority carrier 
lifetime in the material is not yet determined. 
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CONCLUSIONS 


It is difficult to compare the computer calculations with experimental 
results at this point since neither aspect of the research is completed. In 
addition, our measurements have been made on material with large carbon and 
oxygen concentrations and indeed the carbon interstitial and the carbon- 
oxygen-vacancy complex have been the dominant defects observed. The aluminum 
defects predicted by the computer calculations have been observed in EPR 
studies and we are continuing to study the as yet unidentified DLTS peaks in 
aluminum-doped silicon. We are also studying p-type silicon with low carbon 
and oxygen concentrations and with dopants other than aluminum or boron. 

These experiments should be useful for comparison with the damage production 
as modeled by the computer calculations and will also provide information on 
defect characteristics and identity necessary for the construction of the 
model. This joint approach to the problem of radiation damage in silicon 
solar cells should give us the necessary information to produce radiation 
resistant devices. 
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TABLE I.- DEFECTS OBSERVED BY DLTS IN BORON-DOPED 

Peak Energy Level Capture Cross 

Section (cm 2 ) 
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E 

V 

4 
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eV 
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eV 
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eV 
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1 
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0.22 

eV 
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H 9 
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V 

4 

0.20 

eV 



E 1 

E 

c 

- 

0.29 

eV 



A&-DOPED ELECTRON IRRADIATED SILICON 


Annealing Identity 

Temperature 

out 270°C [V 2 ] + 

out 30°C [C ] + 

in 30°C, out 400°C [V+04€] + 

out 400° C ? 

out 170°C [B^+Oj] (tentative) 

in 170°C, out 400°C [B+O+V] (tentative) 

A£ associated 
7 

in 145°C, out 250°C ? 

in at 400°C ? 

in at 445°C ? 


? 



TABLE II.- RATE EQUATIONS 


at = G " 4lr(D v +D i )R vi ™i] - ^D V [V](2 R VV [V] + R VBi tB I i + R VBs [B s 3) + 


Hii - G + 2G - 4 «VVV V1[I1 - 4 V 1 ’ ®IB ‘V + K IVB '™S ] + 

Z, D b 

R V 2 I [V 2 ]) 


IT" ■ G V 2 - + S.D^IV] 2 


a[B J 

3t =- 4irD I R IB [I][Bg] + 4ir D v R Y B [V] [Bj] - 4irD B ^ B t B xHBg] + 
o X X X o 

4 ’ rD B T R V,B T [B I HV 2 1 " W Ab' ,][, S ] + 4 * D I E IVB c [I1[VB S 1 

1 4» 1 U O 


a[B ] 

~~3t = 4 ^i R ib “ ^V^B ^ B i^ V ^ “ 4irD B T R B T B *- B I-**- B S- 

X X X o 


4 ’ iDb i\ 2 Bi IBiIIV 2 1 


3 [B B ] 

3t ‘'VVb'W 

X X D 


3tVB ] 

“Bi— " “AbI'H's 1 - 4 * D I R IVB ^liTOsl 

s s 


[ ] indicate concentration 
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TABLE III.- NUMERICAL VALUES OF PARAMETERS USED FOR COMPUTER CALCULATIONS 
Capture Radii 

= R w = R ve >1 = = ^Bg = ^I 

' X = Rivb s ' Vs ' 5 A 

Generation rates 

12 2 

G = 5 x 10 /cm /sec generation rate of V,I 
9 2 

G = 5 x 10 /cm /sec 
V 2 

Diffusion coefficients 

-3 2 

= 3.16 x 10 cm /sec 
-7 2 

= 4.15 x 10 cm /sec 

— o 

D = 1.00 x 10 cm /sec 

b i 

D A tl " ° 

Rotating satellite 

— 

flux on for 10 sec (g,G = 0) 

V 2 

flux off for 10 ^sec (G,G = 0) 

* o 
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